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1 INTRODUCTION 



The picture of galaxy formation as a hierarchical clustering 
process (e.g. White & Rees 1978; Fall & Efstathiou 1980) 
has gained considerable observational support over the last 
decade (e.g. Giavalisco, Steidel, & Macchetto 1996; Steidel 
et al. 1998; Abraham et al. 1999). It is thought that evolving 
dark matter haloes set the stage for the subsequent gas cool- 
ing and star formation that ultimately result in the plethora 
of galaxies we see in the present epoch. However, the lack 
of a detailed physical theory of how stars form in early, 
evolving environments remains a pressing problem. For ex- 
ample, when gas dynamics and star formation are ignored, 
distributions of hundreds of dark matter satellites around 
larger galaxies have been produced in numerical simulations 
(Klypin et al. 1999; Moore et al. 1999) whereas only tens 
of satellites are observed. In addition, simulations of galaxy 
formation in the context of Cold Dark Matter (CDM) dom- 
inated universes assume that thermal cooling-function pre- 
scriptions can be used to predict how and when gas converts 
to stars. If it is supposed that the bulk of the cooled gas 
goes directly into star formation, then the inescapable high 
cooling rate that occurs in gas within the growing dark mat- 
ter haloes when the universe was much smaller and denser 
leads directly to the so-called 'star formation catastrophe'. 
Far too high a proportion of baryons end up in stars in such 
simulations compared to the observations. Simulations also 
repeatedly fail to produce Milky Way types of galaxies be- 



ABSTRACT 

In an exploration of the birthplaces of globular clusters, we present a careful exam- 
ination of the formation of self-gravitating gas clouds within assembling dark matter 
haloes in a hierarchical cosmological model. Our high-resolution smoothed particle 
hydrodynamical simulations are designed to determine whether or not hypothesized 
supergiant molecular clouds (SGMCs) form and, if they do, to determine their phys- 
ical properties and mass spectra. It was suggested in earlier work that clouds with a 
median mass of several 10*Mq are expected to assemble during the formation of a 
galaxy, and that globular clusters form within these SGMCs. Our simulations show 
that clouds with the predicted properties are indeed produced as smaller clouds collide 
and agglomerate within the merging dark matter haloes of our cosmological model. 
We find that the mass spectrum of these clouds obeys the same power-law form, 
dN/dM (X M~^-^^°-^, observed for globular clusters, molecular clouds, and their in- 
ternal clumps in galaxies, and predicted for the supergiant clouds in which globular 
clusters may form. We follow the evolution and physical properties of gas clouds within 
small dark matter haloes up to z = 1, after which prolific star formation is expected 
to occur. Finally, we discuss how our results may lead to more physically motivated 
"rules" for star formation in cosmological simulations of galaxy formation. 

cause their gas condenses too quickly to be torqued up to 
the observed angular momentum of gas in spirals (Navarro 
& Steinmetz 1997; Weil, Eke, & Efstathiou 1998; Eke, Efs- 
tathiou, & Wright 2000, hereafter EEW). It is hoped that 
the remedy for these ills lies in the addition of the feedback 
of star formation upon the gas. Vigorous supernova activity 
may heat up the gas enough for it to remain fairly diffuse 
and extended (e.g. McKee & Ostriker 1977, McKee 1995). 



The purpose of the present paper is to apply the in- 
sight of star formation processes in the current epoch to 
help solve this problem. The observations of star forma- 
tion in nearby molecular clouds and nearby extragalactic 
systems clearly show that the prescriptions for star forma- 
tion imposed in the cosmological models are too simple. 
Firstly, self-gravitating (molecular) clouds are not supported 
by thermal pressure, but rather by a combination of tur- 
bulent and magnetic pressure (McKee et al. 1993; Pudritz 
2000). Thus, simple thermal cooling time arguments are ir- 
relevant for predicting how molecular clouds convert some of 
their mass into stars. Secondly, stars form in clusters within 
clumpy subregions of molecular clouds, and the total mass 
of these embedded clumps is only a few percent of the mass 
of the entire cloud. Hence, star formation proceeds only in 
sub regions within molecular clouds that are sufficiently self- 
gravitating, or equivalently, are at sufficiently high pressure. 
Thirdly, the mass spectra of molecular clouds and that of 
their clumps have been found to obey an approximate power 
law scaling 
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dN/dM oc M-i-''±°-i (1) 

(eg. Blitz 1993, Williams et. al. 2000). The fact that 
stars form within such special sub-structures in the self- 
gravitating clouds almost certainly has important impli- 
cations for the formulation of more physically motivated 
"rules" for star formation. 

Do characteristics of galactic star formation in the 
present epoch carry over to the epoch of galaxy assembly? 
Strong indications that the answer is yes come from the old 
halo globular cluster systems that have been studied around 
more than 100 galaxies (see review; Harris 2000). Harris and 
Pudritz (1994, hereafter HP) showed that the mass spectrum 
characterizing globular cluster systems around a variety of 
galaxies has nearly the identical form as the mass spectrum 
for local giant molecular clouds (GMCs) and their clumps. 
Because the efficiency of conversion of clump gas to stars in 
bound clusters must be high, the mass spectrum of bound 
clusters should reflect that of their progenitor clump spec- 
trum. Observations show that globular clusters more mas- 
sive than 10^ Mq and more than a few kpc from the cen- 
ters of galaxies survive evolution in their host galaxies rela- 
tively intact (e.g. Harris, Harris, & McLaughlin 1999). This 
is born out by semi-analytical studies of the tidal evapora- 
tion and disruption of globular clusters in isothermal dark 
matter haloes for a very wide range of cluster masses, con- 
centrations, and orbital parameters (e.g. Capriotti & Hawley 
1996). Such calculations show that clusters in the mass range 
lO'' to IO'^Mq can survive for a Hubble time because they 
are massive enough to avoid evaporation but not too mas- 
sive to undergo rapid orbital decay though dynamical fric- 
tion. Thus, the mass spectrum of the more massive globular 
clusters around any galaxy essentially preserves the mass 
spectrum of its progenitor clumps. The numerical model of 
McLaughlin and Pudritz (1996, hereafter MP) showed that 
such clump mass spectra arise naturally through an agglom- 
eration/fragmentation process, wherein clumps and clouds 
build up by the collision and merger of smaller ones within 
any given dark matter halo, and are destroyed as a conse- 
quence of star formation within them. 

Globular clusters also provide empirical guidance as to 
how stars formed during galaxy assembly in a wide variety 
of systems and environments. The globular cluster luminos- 
ity function or, equivalently, mass spectrum, around galax- 
ies has been measured in a variety of galaxy types such as 
dwarf ellipticals (e.g. Durrell et al. 1996), as well as spirals 
and ellipticals. The mass spectra are found to be identical to 
equation (1) within the errors. The mass spectra of globular 
cluster systems is also independent of the mctallicity of the 
host galaxy (e.g. MP). Moreover, globular clusters arc ob- 
served to form today in galaxies that have abundant supplies 
of gas (e.g. Holtzman et al. 1992; Whitmore & Schweizer 
1995; Zepf et al. 1995; Schweizer et al. 1996; Whitmore et 
al. 1999). Thus, globular cluster formation has a universal 
character which requires only a sufficient supply of gas, and 
is independent of host galaxy conditions and details about 
gas metallicity. 

Previous work on this model for cluster formation fo- 
cussed on processes occurring within some given dark mat- 
ter halo, and did not fully embed these ideas within a full 
cosmological framework. Thus, HP and MP constructed 
a model in which globular clusters form within massive 



clumps of 10^ - 10'' M0 that arc embedded within their 
yet larger host, supergiant self-gravitating (partially molecu- 
lar) clouds, denoted SGMCs. The hypothetical SGMC mass 
range - lO'' - W^Mq - is in fact observed for molecular 
gas within merging or starburst galaxies (e.g. Wilson et al. 
2000). In the model, the mass spectrum of the SGMCs (and 
their substructure - the clumps ) arises from the agglom- 
eration of gas clouds by cloud-cloud collisions within the 
evolving dark matter haloes. The model envisages gas pro- 
cesses within any given, static dark matter halo but does not 
include the mergers of such haloes that arise in any hierar- 
chical model for galaxy formation. Thus, globular cluster 
formation is predicted to occxir much as star cluster forma- 
tion does in molecular clouds at the present day, with the 
only difference being that huge gas reservoirs were available 
(through gas infall or ongoing mergers) at these higher red- 
shifts in comparison with current interstellar media. 

In the present paper, we undertake, for the first time, 
a careful examination of the formation of self-gravitating 
clouds that build up during the assembly dark matter haloes 
in hierarchical cosmological models. Our purpose is to de- 
termine whether or not supergiant clouds form, to measure 
their mass spectra, and to determine how they evolve from 
early redshifts to z = 1. Our high resolution simulations 
focus, in particular, on the study of these progenitor su- 
pergiant clouds before the "star formation catastrophe" oc- 
curs (typically at redshifts z ~ 1). We specifically designed 
our simulations to have the resolution needed to study the 
formation of ensembles of individual, self-gravitating (and 
partially molecular) clouds that are the actual seats of star 
formation. The idea is to test and survey the initial con- 
ditions for star formation, before problems such as the star 
formation catastrophe become endemic. This procedure will 
allow us to test current models of star and globular cluster 
formation on cosmological scales. It will also ultimately lead 
to the development of more physically motivated rules for 
stellar formation in cosmology. 

Our results clearly show that SGMCs are commonly 
formed within evolving dark matter haloes in our first CDM 
cosmological model, and that they have the the mass spectra 
and cloud properties predicted for clouds capable of form- 
ing globular clusters. We discuss our numerical methods in 
§2, and present our results in §3. Implications of our work 
for star formation during galaxy assembly and, hence, for 
globular cluster formation follow in §4. 



2 NUMERICAL METHODS 

We investigate the formation and evolution of star forming 
clouds in one of the several viable, spatially flat, CDM cos- 
mological models (Thomas et al. 1998), known as tCDM. 
Structure formation in dark matter and gas is followed us- 
ing smoothed particle hydrodynamical simulations (SPH). 
We evaluate the formation of supergiant star-forming clouds 
during the evolution of small density perturbations initially 
seeded in an expanding universe. 

In order to attain the mass resolution required to iden- 
tify massive, self-gravitating clouds and their substructure, 
we choose small dark matter haloes from a large cosmolog- 
ical simulation which has been evolved to the present day 
{z = 0). After tracing the dark matter haloes back to red- 
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shift z = 24, their resolution is improved by increasing the 
particle numbers by factors of hundreds and adding an equal 
number of gas particles. The volume outside the high res- 
olution region is represented by high-mass, low resolution 
particles to retain the full gravitational characteristics of the 
original model. The models are then evolved from z = 24 to 
z = 0. Section 2.1 presents properties of the large cosmolog- 
ical simulation and the method of producing high resolution 
models. The properties of the dark matter plus gas simulsr 
tions the code used to evolve the high resolution models are 
presented in §2.2. 

2.1 Cosmological Model 

A dark matter only cosmological simulation is used to gener- 
ate initial conditions for the high-resolution dark matter plus 
gas simulations described below. The cosmological model is 

one of two spatially flat CDM simulations by EEW. 

The tCDM cosmology used in this paper has a mat- 
ter density parameter Qrn = 1, zero cosmological constant, 
and a Hubble constant h = 0.65 in units of 100kms~^Mpc~'^. 
The power spectrum that generates the scale-free initial den- 
sity fluctuations for these simulations has a shape parame- 
ter of r = 0.2, which fits observed large-scale correlations 
in galaxy and cluster surveys (Efstathiou, Bond, & White 
1992). This specific cosmological model could bo relevant to 
the situation in which the late decay of r neutrinos produces 
an additional background of relativistic e~ and n neutrinos. 
This has the effect of delaying the onset of matter domi- 
nation. The other large cosmological simulation produced 
by EEW is a ACDM model with n„, = 0.3 and Qa = 0.7. 
The ages for the rCDM and ACDM cosmological models arc 
10.0 Gyr and 14.5 Gyr, respectively. Due to its shorter evolu- 
tion time, we chose to explore the formation of star-forming 
clouds in the rCDM first. No significant differences in the 
physical processes related to the evolution of SGMCs arc 
expected for other cosmologies, although the precise time 
when SGMCs first turn up will of course be different (see 
§4 for further discussion). The results on the formation of 
star-forming clouds in ACDM models will be reported in a 
subsequent paper. 

The large dark matter simulation represents a cubical 
volume of the universe of length L = 32.5/i~^ Mpc with 
N = 128^ dark matter particles. Figure 1 shows a 2-Mpc 
deep slice of the simulation at z = 0. EEW identified viri- 
alized haloes using a group-finding method, locating over 
2000 haloes with particle numbers ranging between about 
40 and 30,000. The largest dark matter potential wells are 
easily seen at the nodal intersections of the filaments. One 
of the small haloes we chose to resimulate at high resolution 
is identified by the arrow near the top middle of the figure. 
This halo contains 65 dark matter particles and has a mass 
M = 4.5 X IO^^Mq. 

Haloes to be rcsimulated at higher resolution are chosen 
using a set of balanced criteria. Firstly, the mass of each halo 
is required to be relatively small because computing time 
increases as the number of particles increases. We represent 
the high resolution region of the simulation by Ndm = 34^ 
dark matter particles and Ng = 34'^ gas particles. The sur- 
rounding low resolution region contains approximately 5000 
dark matter particles representing the gravitational forces 
due the the mass outside the high resolution region. In or- 



der to achieve the mass resolution necessary to examine the 
structure supergiant star-forming clouds, the masses of the 
small haloes are limited to approximately 5.0 x lO'^'^ Mq. 
Secondly, we choose haloes comprised of enough particles to 
be physically believable structural entities. The large cos- 
mological simulation is searched for haloes containing « 60 
particles. When expanded to high resolution, the haloes con- 
tain gas particles with masses between 2 and 5 x 10® Mq. 
Thus, while we have sufficient numerical resolution to study 
the formation of SGMCs, our simulations do not yet allow 
us to resolve the individual clumps within these clouds that 
will be the actual sites for the formation of star clusters. In 
these simulations, the baryon density parameter Qi, = 0.06; 
thus, 94% of the mass is contained in the dark matter par- 
ticles. Finally, haloes are required to be well-separated from 
any other large mass condensations so that computational 
power will not be exerted upon nearby regions of little in- 
terest. 

The high resolution initial conditions are created using 
the method described by Weil et al. (1998). After a halo 
is chosen at 2 = 0, all particles within a spherical volume 
surrounding it are located at z = 24, the redshift at which 
the cosmological simulations commence. A cube just large 
to contain all the particles is constructed. The length of the 
cube is typically a few Mpc. Populating the high resolution 
cube requires displacing particles from a regular lattice in 
order to represent the spatial fluctuations in density of the 
cosmology. In order to correctly populate their finer vol- 
umes, these high resolution particle displacements require 
additional calculations that add power above the Nyquist 
frequency of the original simulation. The higher mass, low 
resolution region particles are then added in a spherical dis- 
tribution. 

2.2 Numerical Simulations 

In order to allow statistical analysis of the properties and 
evolution of star-forming clouds, seven different dark mat- 
ter haloes were chosen for high resolution resimulation. We 
studied the evolution of each of these systems, which we 
designate by distinct Run nunil)ers: 1 7. Each Run was 
evolved using the TREESPH code (Hernquist & Katz 1989). 
This code treats both dark matter, which is purely gravitar 
tional, and gas, which includes hydrodynamical and radiar- 
tive physics (for details, see Weil et al. 1998). 

The gravitational forces between particles are computed 
using a hierarchical tree algorithm (Bames & Hut 1986). 
Distant particles are grouped into hierarchical cells in which 
the potential is approximated using a truncated multipole 
series. The forces on each particle are updated by examin- 
ing each cell in the tree. The ratio of the size of the cell to 
its distance from the particle is compared to a defined tol- 
erance parameter to determine whether the cell should be 
subdivided or the forces due to it calculated as a whole. In 
this way, computational time is reduced by avoiding explicit 
two-body force calculations for remote particles, while an 
accuracy of ~ 0.1% is maintained by a tolerance parameter 
of 0.6. 

The gravitational potential is softened using a cubic 
spline in order to reduce two-body relaxation. The softening 
lengths assigned to the particles define the length resolution 
of the code. For dark matter particles, the softening edm 
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Figure 1. 2-Mpc deep slice of the large, dark matter only cosmological simulation of EEW. The arrow near the top middle points to 
the small halo which is resimulated at high resolution and shown in Figures 2 and 3. 



is always between 2 and 3 kpc. For the gas partieles, with 
masses more than a deeade smaller than the dark matter, tg 
is between 0.6 and 1.0 kpc. This gaseous softening length im- 
poses another restriction on the size of clouds that we can 
physically resolve, namely, nothing less than several kpc. 
Given that the HP model suggests a fiducial value of one 
kpc for the median SGMC in the distribution, we caution 
that our simulations will inevitably produce fluffier, lower 
density clouds by virtue of the limitations imposed by the 
softening length. Detailed analysis of cloud structure will 
require future, higher resolution simulations. 

Hydrodynamical interactions between gas particles are 
treated using smoothed particle hydrodynamics (Lucy 1977, 
Gingold & Monaghan 1977). The gas is partitioned into 
self-gravitating fluid elements that evolve according to La- 
grangian hydrodynamic conservation laws. Explicitly, pres- 
sure gradients and the viscosity of the medium are calculated 
using conservation of momentum: 



and 

^ = VP + a^'^^-V^W, (3) 

where r and v are vector positions and velocities, pgas is the 
gas density, P is the pressure, a""'^ comprises the acceler- 
ation due to viscosity, and 4> is the gravitational potential. 
The equation of state if that of an ideal gas; stability dur- 
ing coUisional shocks is maintained by an artificial viscosity 
(Navarro & Steinmetz 1997). 

In order to treat the evolution of the thermal energy, 
radiative cooling is computed from the ionization, recom- 
bination, and cooling rates of Black (1981, Table 3) with 
modifications by Cen (1992). Analytic approximations to 
observed rates of collisional ionization, recombination, col- 
lisional excitation, bremmstrahlung, and Compton heating 
and cooling are constructed for H I, H II, He I, He II, and 
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He III. Because we ignore star formation, rather than ini- 
tiahzing and evolving metaUicities for the gas particles, pri- 
mordial gas is assumed to have a hydrogen mass fraction of 
X = 0.76 and a mass helium fraction of y = 0.24. As already 
noted, observations of the IMF and the overall characteris- 
tics of the mass functions of globular clusters indicate they 
are independent of the metallicity of the gas out of which 
they formed. 

The local quantities appearing in the hydrodynamical 
equations are computed by interpolating between gas parti- 
cles (Monaghan 1985, Steinmetz & Miiller 1993). Individ- 
ual smoothing lengths are calculated for each particle in 
order to ensure that the values for each particle arc cal- 
culated with reference to 25 to 45 neighbors. In addition 
to individual smoothing lengths, individual timesteps are 
permitted for each particle. The largest time step for these 
TREESPH simulations is At = 2.5 x 10^ years; the smallest 
number of timesteps in which a particle can reach z = is 
Nstep = 4000. Owing to the need to satisfy local stability as 
required by the Courant condition, particles in the highest 
density areas are allowed to reduce the size of their timesteps 
by up to a factor of 32. In simulations with a large range 
of mass inhomogeneitics, this adaptivity in space and time 
permits areas of high density to be resolved to an arbitrary 
resolution. 

Our simulations deliberately avoid the implementation 
of any star formation rule. This enables us to study the 
properties of ensembles of self-gravitating clouds that build 
up before any star formation catastrophe - predicted by a 
naive interpretation of gaseous thermal cooling processes - 
could be important. This assumption implies that we can- 
not chart the evolution of the metallicity of these early sys- 
tems before z = 1. We therefore refrain from speculating 
about the observed metaUicities of globular cluster systems, 
or any observed trends of this kind (e.g. age-metallicity re- 
lations) in any given galaxy, and how those may come to bo. 
Without the high spatial and temporal resolution necessary 
to resolve and follow the star formation in clumps within 
SGMCs, such discussion would be premature. The second, 
and related, limitation is that the destruction of clouds due 
to star formation is not taken into account. The proper treat- 
ment of cloud destruction critically depends however, upon 
the incorporation of correct star formation physics. Both of 
these limitations will be relaxed in subsequent work. 

The quantities discussed thus far refer to the fiducial 
initial conditions for our simulations. However, each of the 
seven haloes was first resimulated at "low resolution." The 
low resolution simulations were assigned particle numbers 
for dark matter and for gas of Niow ~ 22'^ each in compar- 
ison to the Nhigh = 34"' each of the high resolution simu- 
lations. Simulations at this resolution require less than one 
day to compute on a Sun Ultra 10 workstation. Each halo 
was evolved from a = 24 to z = in low resolution; and was 
examined for any aberrations that would bias the results. 
Systems that were in the process of merging or consisted of 
several subsystems at a = were excluded. Subsequently, 
each halo was evolved in high resolution from 2; = 24 to 
z = l. 

Figures 2 and 3 show the evolution of dark matter and 
gas, respectively, during Run 2. The length of each panel 
is 2 Mpc in comoving coordinates and the redshifts are 
marked at the upper right. In this hierarchical cosmology. 



several dark matter potential wells, into which star-forming 
gas falls, have formed by z = 5. The gas is more centrally 
concentrated than the dark matter, forming small, denser 
clumps, surrounded by more diffuse gas. If the first epoch 
of efficient star formation occurs at z = 5 in our adopted 
tCDM cosmology, then the present age of the first objects 
is 9.7 Gyrs. In comparison, if stars first form efficiently at 
this time in a standard CDM cosmology, their age is w 11 
Gyrs. Even at « = 1, several gas condensations, which will 
eventually merge to form a small galaxy, are visible. The 
present age of stars formed at this time is 6.5 Gyrs in tCDM 
and w 8 Gyrs in standard CDM. We discuss the significance 
of these ages in §4. 



3 RESULTS 

In order to construct a detailed theory for the formation 
of globular cluster systems (GCS) and, more generally, for 
star formation within hierarchical cosmologies, the results of 
simulations need to be compared to observational quantities 
and models. If stars do indeed form at early epochs in the 
high pressure clumps within massive gas clouds (as they do 
in the present epoch) , then high resolution cosmological sim- 
ulations of galaxies should enable the identification of these 
supergiant clouds. Herein, bound gas masses that represent 
the star-forming clouds within Runs 1-7 are found and ex- 
amined. The time evolution and statistical and individual 
properties of these supergiant, star- forming clouds are also 
investigated. 

Our main goal was to identify sufficiently massive bound 
gas clouds formed in high resolution simulations such as 
that shown in Figures 2 and 3. Two methods of identify- 
ing clumps of particles were applied to outputs at several 
redshifts of Runs 1-7. Because the clumps interact with 
one another and can be irregular, overlapping, and not in 
dynamical equilibrium, we require an algorithm which can 
uniquely locate separate clumps and ensure that the systems 
are bound. 

The spherical overdensity algorithm (SO; Lacey & Cole 
1994) grows spheres around density maxima until they 
reached a defined density contrgist with the background. The 
overdensity A = 178 corresponds to virialized regions within 
a uniform collapsing sphere. The friends-of-friends algorithm 
(FOF; Davis et al. 1985) uses a linking method by which 
particles are identified as members of a clump if they are 
closer to the nearest member of the clump than a defined 
'link length.' In essence, the link length determines a local 
density criterium. The results of the two clump-identifying 
algorithms supports the Lacey and Cole (1994) observation 
that small clumps are sometimes merged with larger ones for 
the fiducial parameterizations. The fiducial FOF link length 
is 0.2 times the mean interparticle separation. Herein, how- 
ever, the typical link length is defined to be about a factor 
of two less than that, in order to separately identify small 
clumps and to avoid including unbound particles. It is likely 
that, for many clumps, several outlying particles are not in- 
cluded Eis members; however, the mass contribution of these 
few particles is negligible in most cases. 

For the several output redshifts of Runs 1- 7, all clumps 
with particle numbers A^c.g > 15 and Nc,dm > 30 for gas and 
dark matter, respectively, are located. Subsequently, only 
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Figure 2. Time evolution of Run 2 dark matter. Each panel measures 2 Mpc in length in comoving coordinates. Redshift is shown in 
the upper right hand corner of each panel. 



those clumps whose particles arc enclosed within the virial 
radius of the small galaxy at z=0 arc retained. Those clumps 
with gas masses greater than a few times IO^Mq are con- 
sidered to be supergiant star-forming clouds. We apply the 
clump-finding algorithm to three different types of bound 
objects, namely, to gas plus dark matter agglomerates, to 
pure gas clouds and, finally, to pure dark matter haloes. 



3.1 Cloud Populations 

The epochs at which star-forming clouds appear are inves- 
tigated for Runs 1 - 7. At our mass resolution limit of a 
few times 10® Mq, by z = 5 each of Runs 1-7 produces 
2 to 10 gas clouds with masses of several times 10^ to sev- 
eral times 10**A/q. These will eventually become part of the 
small gala^xy. Figure 4 shows the evolution of supergiant star- 
forming clouds located by FOF at eight redshifts for Run 2. 
Each of these gas clouds consists of particles identified to 



be contained within the virial radius of the small galaxy at 
z = 0. Each panel measures 1 Mpc in length. In the bottom 
right panel, the gas particles at z = are shown. The circle 
outlines the axea within the virial radius of the galaxy, which 
encloses 4.4 x IO^^Mq (including dark matter and gas). At 
2 = 1, the clouds contain 1.2 x 10^°Mq gas. 

Figure 4 shows the bound gas clouds at eight differ- 
ent redshifts, and their final product at 2 = 0. The number 
of star-forming clouds increases with redshift eis the den- 
sity perturbations of the hierarchical cosmology grow large 
enough to contain bound material. Although initially well- 
separated in space, clouds flow along the filaments visible in 
Figures 2 and 3 towards the dominant density concentration, 
piling up in the center as they merge to form a galaxy. The 
clouds appear irregular in shape throughout all redshifts. It 
is clear from this figure that clouds are growing by cloud- 
cloud collisions at the same time as the galaxy is building 
up through a hierarchical merger process. 
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Figure 3. Same as Figure 2 except gas is shown. 



One can clearly see the collision and agglomeration of 
clouds in Figure 4. At z = 3.0, wo sec 3 distinct clouds in 
middle of the left side that collide, leaving two clouds at 
that position in the panel at z = 2.5. As another example, 
the two amorphous clouds in the upper middle paxt of the 
z = 2.5 panel collide to produce the larger cloud at the same 
position in the z — 2.0 panel. Finally, the group of clouds in 
the middle of the panel a.t z = 2.5 collides and agglomerates 
hy z = 1.5. These, and other details clearly show that gas 
agglomerates into more massive structures by cloud-cloud 
collisions. Gets is brought together in two ways: 1) it falls 
into an individual dark halo and collides with other clouds 
in that halo, and 2) gas in a given halo collides with the 
fresh gas contained in other small dark matter haloes that 
are merging. 

Figure 5 shows the central region of Run 2 at 2: = 1. 
The left panel shows the gas and the right panel shows the 
bound gas clouds. The colors of the clouds are non-physical; 
all the gas particles that belong to a distinct cloud are rep- 



resented by a single color. The dark matter particles are not 
included in this view. As a rule, clouds consist of a dense 
central core, which contains most of the mass, surrounded 
by more diffuse gas. We caution that we cannot further re- 
solve the expected substructure within our clouds, because 
of the limitations imposed by our resolution. In the same 
way that the halo in the original cosmological simulation of 
Figure 1 is essentially featureless, these clouds, in which the 
mass of each gas particle is on the order of the mass of a 
globular cluster, are not complex. However, if we examine 
the ensemble of clouds in the center of Figure 5 as a whole, it 
has the appearance of an amorphous agglomerate containing 
several high density cores. If the mass resolution for Runs 1 
- 7 is further increased, the individual clouds would likely 
show further geometrical complexity. 

Each of the seven separate evolving haloes is accompa- 
nied by a spectrum of masses for both the gas clouds and 
the dark matter. Because the number of clouds in each sep- 
arate Run can be small, the collections of bound objects in 
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Figure 4. Bound gas clouds at eight redshifts between z = 5 and z = 1 for Run 2. The clouds merge eventually to become part of the 
small galaxy shown at z = 0. The circle outlines the area within the virial radius. Each panel measures 1 Mpc in length in comoving 
coordinates. 



each of Runs 1-7 are added together at each redshift in 
order to give a larger statistical sampling. At a redshift of 
z = 5, the haloes consist of between 2 and 11 clouds each. 
At intermediate redshifts, most haloes consist of about 15 
clouds but one contains only 3 and another contains about 
30. At 2 = 1, half of the haloes contain about 6 clouds be- 
cause many clouds have merged, one contains about 15, and 
the rest have about 30. At a = 5, the total number of clouds 
used to calculate the mass spectrum is 30; at later redshifts, 
the total number is greater than 100. Masses are calculated 
for each cloud. Figure 6 shows mass spectra for all the gas 
clouds at nine redshifts. Our data appear to flatten slightly 
at the low mass end, which may be a consequence of our 
resolution. Remarkably, least squares fits to the logarithmic 
data show that, at all redshifts, a power law fits the data very 
well for masses > 10^'^Mq, with 1.6 < a < 1.9. There is, 
at best, only very weak evidence for an evolutionary trend 
for a decreasing over time. The figure indicates that the 



original low mass clouds accrete mass or merge with other 
clouds, moving to the right in the panels over time. The low 
mass clouds are continually replenished as new density con- 
densations grow. The low mass end of the spectrum is not 
depleted but contains essentially the same number of clouds 
at all redshifts. Thus, the power law slope is expected to be 
steeper at high redshift, before high mass clouds have had 
time to build up, but flattens over time as mass moves to the 
right in the figure. Because clouds are simply agglomerating 
without being shredded by feedback from internal stax for- 
mation (the latter effect was included in MP's model), the 
upper mass end of the spectrum coritirmes to be pushed to 
higher masses. It is likely that a correct treatment of star 
formation within these clouds would result in the truncation 
of this mass spectrum. At this stage, therefore, the predic- 
tions of the analytic HP model and its numerical extension 
(MP) are seen to arise very naturally within a hierarchical 
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Figure 5. Run 2 at z = 1. The left panel shows the gas in the central region and the right panel shows the bound gas clouds. Each 
panel measures 0.5 Mpc in length in comoving coordinates, as labeled on the left and bottom axes. Also shown are physical Mpc on the 
right and top axes. 



cosmological model. We discuss the ramifications of this in 
§4. 

Figure 7 shows the mass spectra for the dark matter 
in the bound objects. The solid line is the least squares fit. 
Again, the GCS power law fits well, with 1.5 < Q < 1.8. The 
evolutionary trend of the slope is weaker here. If the trend is 
cosmological, one explanation is suggested by Figures 2 and 
3. The dark matter forms potential wells into which gas falls, 
as seen by comparing the number and density of clumps at 
any one early redshift, say z = 3, in the two figures. Once 
the dark matter clumps are formed, they do not increase 
in mass as the associated gas does. Figure 7 shows that 
the mass does not evolve strongly to the right over time, 
but that most high mass dark matter clumps are in place 
hy z — 3. The fiattening or turnover at small mass is even 
more apparent for the dark matter than the gas spectra. This 
turnover also occurs in the globular cluster systems of real 
galaxies (MP). In observed GCS, it is difficult to determine 



whether the change in the power law slope is physical or 
due to incomplete sampling at low masses. In Figure 7, it 
appears that a collection of small dark matter clumps, with 
masses less than Mq, form at intermediate redshifts but 
higher resolution simulations are needed to determine if the 
turnover is physical or due to limited resolution. Meanwhile, 
we limit the fit to the mass spectra to clumps with masses 
greater than 3 x 10** Mq. 

The dark matter mass spectra in Figure 7 are also com- 
pared to the Press-Schechter linear analytic description of 
the evolution of mass in a hierarchical cosmology (Press & 
Schechter 1974). The Press-Schechter multiplicity function 
predicts how bound objects condense out of a density field 
described by the cosmological power spectrum 15^1 cx fc", 



N{M, z)dM oc M 



n/6-3/2 



exp 



M 



dM (4) 
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Figure 6. Mass spectrum of combined gas clouds for all seven runs shown at nine redshifts. The solid lines show the power-law fits of 
eq. (1). The slopes, a = Alog{dN / dM) / AlogM , are shown at the top left, below the redshift. 



where M,{z) oc (1 + z)* '*''("+-^)\ The power spectrum of 
our cosrriological simulations is a scale invariant spectrum 
with n — 1 but extrapolated to smaller scales (Efstathiou, 
Bond, & White 1992). Figure 7 shows fits using a non- linear 
Marquardt method to masses M > 3.5 x 10^ Mq. The dotted 
and dashed fines are for n = 1 and n = 0, respectively. 

The Press-Schechter method is limited by the assump- 
tion of linear growth of bound objects but can be extended to 
treat non-linear epochs of collapse and merging (e.g. Bower 
1991, Bond et al. 1991) Application of the original Press- 
Schechter mass function to our data shows that the fit di- 
verges at high masses where the exponential part of equa- 
tion (4) is important. In addition, the reliance on redshift in 
the exponential is weak in the data. Below the exponential 
cutoff, the fit is a simple power law. However, for n = 1, 
the expected slope is —1.3, fiatter than the least squares 
fit slopes. That agglomeration-produced, power-law models 
fit the dark matter haloes better than Press-Schechter indi- 



cates that the effects of merging overwhelm those of linear 
growth, as expected in hierarchical cosmologies. 



3.2 Individual Cloud Properties 

Whereas the mass spectrum predicted by observations of 
globular cluster systems is fit at all redshifts by the super- 
giant star-forming clouds assumed to be their birthplaces, 
individual properties of clouds are partially determined by 
the resolution of the simulations. The HP SGMC properties 
are based on the application of the virial theorem. They pre- 
dicted that median SGMC in the power-law spectrum would 
have a mass ~ 7 x 10* Mq, a radius r ~ 900pc, surface den- 
sity E ~ 260Mq/pc^, and density p ~ O.22M0/pc^. Figure 
6 shows that the low mass end fit by the power law is pop- 
ulated by clouds of the expected mass, a few times 10* to a 
few times 10®. 

The assumptions of approximate virial equilibrium for 
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Figure 7. Mass spectrum of combined dark matter clumps for all seven runs sliown at nine redsliifts. Tlie solid lines show power-law 
fits, as described in the previous figure caption. The dotted and dashed lines are fits of the Press-Schechter multiplicity function of eq. 
(1) with n = 1 and n = 0, respectively. 



the SGMCs can be tested by noting that the virial theorem 
for uniform spherical sclf-gravitatiiig objects predicts that 
a cloud's mass is related to its radius and internal (1-D) 
velocity dispersion by M = 5cr^r/G. The state of our clouds 
is tested by plotting the virial theorem in the form M versus 
the quantity a^r for the numerical data. Figure 8 shows the 
logarithmic expression of the quantities for the clouds of 
Runs 1 -7 at nine redshifts. Least squares fits give slopes, a, 
close to unity, indicating that, as a whole, cloud ensembles 
arc indeed near virial equilibrium. However, Figure 8 also 
shows tfiat some individual clouds are not virialized but are 
just condensing out of the primordial density perturbation 
field or affected by recent interactions (objects of higher than 
average a^r with masses around 10* Mq as an example). 

With gas particle masses of a few times 10^ Mq, the 
length resolution in the simulations is limited to « Ikpc. Be- 
cause the gravitational potential in the simulations is soft- 
ened in order to avoid two-body effects, the radii of the 



clouds do not evolve with their mass. Over three decades in 
mass and between rcdsiiifts z = 5 and z = 1, cloud radii 
range from 0.5 < r < 5 kpc. Due to resolution limits, a typ- 
ical cloud is approximately twice the predicted size. Thus, 
the typical surface density it a factor of four too low and 
the typical density a factor of nine too low to be compared 
directly to SGMCs. 

Any model of hierarchical structure formation that con- 
tains some prescription for cooling the gas produces an inter- 
stellar medium that has both cooler, self-gravitating clouds 
and a hotter, diffuse medium. The diffuse ISM provides the 
pressure that truncates the sizes of molecular clouds in our 
own and other galaxies. Similarly, SGMCs are assumed to 
be isothermal clouds confined by a surface pressure due to 
a hotter external medium (e.g. HP). 

The temperature structure of the star-forming clouds in 
Runs 1-7 cannot explicitly represent that of the interstellar 
medium because star formation, feedback, and other meth- 
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Figure 8. Virial theorem in form M oc a^r for Runs 1 
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7 gas clouds. Redshifts and least squares fits are shown at the top right in each 



ods of heating the ISM are not treated in our simulations. 
Nevertheless, the gas has a temperature structure &t z = 1 
in which the temperature of particles in the centers of the 
clouds is < 10, OOOK, the immediately surrounding material 
is between 10, OOOK and 12, OOOK, around which particles 
of > 12, OOOK are smoothly distributed. As expected, the 
centers of the clouds, where star formation should occur in 
cool, dense cores, contain low temperature gas, while the 
more diffuse surrounding medium is hotter. Although not 
fully representative of the multiphase ISM, the temperature 
structure is that of a clumpy, warm, neutral medium sur- 
rounded by a more diffuse, warm, ionized medium. 



4 DISCUSSION AND SUMMARY 

The simulations clearly establish that, in tCDM cosmolog- 
ical models at least, the hierarchical growth of dark matter 
haloes is accompanied by the growth of self-gravitating, su- 



pergiant clouds through the related processes of cloud-cloud 
collisions and agglomeration. The mass spectrum that we 
find for SGMCs takes the expected power-law form. That 
clouds are more concentrated than the dark haloes indi- 
cates that ISM properties are the key to cloud production. 
Our SGMCs are not merely condensations in the centers of 
mini-haloes. Rather, they are self-gravitating entities that 
are similar to very massive GMCs that form in starburst 
galaxies today. 

We have good indications that these results pertain to 
other cosmological models as well. We evolved one halo from 
the ACDM model mentioned in §2, and found similar mass 
spectra and other features of the SGMCs reported here. The 
generality of our results is, in fact, to be expected on physical 
grounds because, regardless of the hierarchical model for the 
dark matter, the gas dynamical processes remain the same. 
Gas will collapse into the haloes, agglomerate there, and 
agglomerate with the gas in other nearby haloes that merge 
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with their hosts. Gas agglomeration is a robust process that 
has been studied for decades in many ISM contexts (e.g. 
reviews by Elmegreen 1993, Pudritz 2000), beginning with 
Oort's famous model for gas clouds in the ISM. In a future 
paper, results of a study of gas physics and SGMC formation 
in ACDM models will be reported. 

It has often been suggested that star formation within 
assembling galaxies is restricted to the disks of spiral-like 
galaxies. Thus, it is interesting to observe that the spatial 
distributions of SGMCs are not disk-like in any of Runs 1- 7, 
but rather are spread out in an irregular fashion throughout 
the still-assembling dark matter haloes that, at z = 0, will 
bo the haloes of galaxies with masses of a few xlO" M0. If 
SGMCs give rise to star and globular clusters, the clusters 
are likely to be distributed throughout the halo of the galaxy 
at our epoch. 

Several models assume that globular clusters form by 
fragmenting out of massive, self-gravitating sheets of gas. 
These sheets could result from cloud-cloud collisions (e.g. 
Kumai et al. 1993) or swept-up supershells from earlier 
epochs of primordial star formation (e.g. Brown et al. 1995). 
Collisions between clouds, it has been suggested, could com- 
press the gas into thin, dense, post-shock sheets that may 
only be 10s of parsecs thick. Gravitational fragmentation 
is then invoked to form clusters by gravitational instability 
in such self-gravitating layers. Our simulations do not sug- 
gest that our SGMCs are significantly crushed by collisions 
although we hasten to emphasize that our gas smoothing pa- 
rameter would not allow us to resolve any layer thinner than 
0.6 kpc. The central point remains: we find that SGMCs are 
assembled in lower mass potential wells systematically with 
time, as Figure 6 shows. Cloud-cloud collisions in lower mass 
dark haloes favors their agglomeration, rather than their dis- 
ruption, because the cloud-cloud velocity dispersion in lower 
mass haloes is significantly lower than it is in more massive 
galaxies that emerge later. Secondly, we note that a major 
prediction of sheet fragmentation is that there is a preferred 
mass scale that corresponds to the most rapidly growing, 
gravitationally unstable mode in the sheet (e.g. Elmegreen 
& Elmcgrccu 1979). The expected, strongly peaked mass 
spectrum that results from this process is at odds with the 
observed, power-law mass spectrum of clumps and globular 
clusters. We suggest that sufficient mass, which can be as- 
sembled by a variety of different physical processes - and not 
external triggers - is all that may be required to form clus- 
ters in the ISM of the Milky Way, as well as in assembling 
galaxies at high redshift (see e.g. HP). 

Observations of star formation in current GMCs indi- 
cate that the typical star forms as a member of a star clus- 
ter. Most star clusters formed in the molecular clouds of 
our galaxy today become unbound as the natal clumps in 
their host GMCs are dispersed. The same process may oc- 
cur for cluster formation in SGMCs. While we anticipate 
that SGMC internal clumps will be the sites for all of the 
star formation in the assembling galactic halo, only a few 
of the putative clusters will be sufficiently bound to survive 
clump dispersal. Bound clusters become globular clusters. 
Those that do not remain bound - because of lower star 
formation efficiency within their natal clumps - will dis- 
perse (as the open clusters in our own Milky Way do today) 
and produce a halo star population. It is likely that these 
clouds also produce the stars that account for the very faint 



(28mag arcsec"^), stellar luminosity observed at large radii 
around many nearby galaxies (Malin & Hadley 1999). Proof 
of these assertions must await further simulations that can 
include reasonable rules for star formation. 

Which aspects of the formation and evolution of star- 
forming clouds are affected by the underlying cosmological 
model? The most obvious is that the cosmological model 
dictates the rate at which the dark matter haloes evolve. 
This, in turn, determines when the gas undergoes substan- 
tial cooling and agglomeration. Thus, the model determines 
the time at which objects such as globular clusters will be 
produced. The measured ages of the globular clusters can 
provide an important constraint on any simulation of star 
formation during galaxy formation. By requiring that cos- 
mological models form SGMCs and, presumably, globular 
clusters within the age constraints provided by observed 
globular clusters, the field of viable theories is further nar- 
rowed. 

Consider the case of our rCDM simulations. The earli- 
est time at which SGMCs in the right mass range begin to 
appear in the simulations is a redshift of order z ~ 8. This 
is an approximate result; yet higher resolution is needed to 
determine what occurs at higher rodshifts. The oldest that 
a globular cluster could be in this cosmology is ~ 10 Gyr if 
globular cluster formation began a,t z = 8. This is at odds 
with the latest calculations of globular cluster ages. If the 
discrepancies that result from distance scale uncertainties 
are considered, the upper limit of Milky Way globular clus- 
ter ages can range anywhere from 12 to 17 Gyr (Grundahl 
et al. 2000). By invoking further constraints, these authors 
suggest an age of 16 Gyr for the oldest cluster. Even with 
this large latitude, rCDM models with representative Hub- 
ble constants cannot account for globular cluster ages. On 
the other hand, ACDM models, with their longer ages due 
to accelerating expansion, likely provide a reasonably com- 
fortable age to admit globular clusters. 

Nevertheless, the lack of feedback from star formation 
on the process of SGMC assembly is an important limita- 
tion. A steep drop in the mass spectrum of globular clus- 
ters above several million solar masses is interpreted by 
HP as arising from the destruction of SGMCs through su- 
pernova activity. A steady-state solution to the agglomera- 
tion/fragmentation equations arises if cloud growth is bal- 
anced by the destruction of the most massive clouds through 
disruption by star formation (eg. Oort 1954, Kwan 1979). 
When the time-scale for forming a yet more massive cloud 
becomes longer than the cloud's destruction time due to 
massive stars and supernovae within it, a strong break in the 
mass spectrum along with a sharp reduction in the number 
of clouds is seen. 

A steady state mass spectrum for the SGMCs does not 
occur in our models because star formation is ignored. Thus, 
our SGMCs establish the correct power-law form, and con- 
tinue to grow to ever larger masses throughout the simula- 
tion. The disruption of SGMCs must occur at the very latest 
by 2 = 1, when all indications are that a great deal of star 
formation must occur. However, globular cluster data sug- 
gest that this truncation occurred at higher redshifts when 
supernovae ripped their natal SGMCs apart (MP). Prescrip- 
tion of a correct rule for star formation will inevitably trun- 
cate the mass spectrum of SGMCs. 

The fundamental limitation of the simulations pre- 
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sented in this paper is that we are unable to resolve the 
substructure of the SGMCs. Star formation will take place 
in the clumps within our clouds, but the expected masses 
of the clumps is less than the minimum mass of our simu- 
lation particle. As wc have made abundantly clear however, 
it is hazardous to blindly apply star formation rules to the 
bulk of the gas in an SGMC. An analysis of the onset of 
star formation must await a future, much higher resolution 
simulation. Nevertheless, we have clearly shown that we can 
resolve and study the SGMCs which provide the arenas in 
which clump, and cluster formation, are likely to take place. 

A limitation to all current simulations is that the sup- 
port of clouds by magnetic pressure, as well as by MHD 
turbulence of some kind, is not modeled in the purely hy- 
drodynamic codes. This is not as irrelevant as might first ap- 
pear. Any local molecular cloud undergoes rapid cooling as 
it forms dues to the efficiency of CO cooling. But, as already 
noted, this does not lead to catastrophic rates of star for- 
mation (i.e. within a dynamical time) throughout the cloud. 
Catastrophic star formation is prevented because of the ef- 
fective support provided by MHD turbulence and magnetic 
pressure within the clouds. Thermal pressure becomes ir- 
relevant. Observations indicate that dynamically significant 
magnetic fields are present in high redshift galaxies (Kro- 
nberg 1994) and, therefore, are an important factor in the 
physics of any self-gravitating clouds found within them. 
Magnetic fields could be generated in the turbulent, collid- 
ing gas in the Eissembling dark haloes by localized dynamo 
processes (e.g. Pudritz & Silk 1989). Magnetic fields of the 
required strength (several jj, Gauss) can predate the forma- 
tion of the first stars and should also be addressed within 
improved numerical models. 

In conclusion, we have found that the supergiant self- 
gravitating clouds (SGMCs) do indeed form in the dark 
matter haloes in our rCDM simulations. The SGMCs as- 
semble through cloud-cloud collisions as gas collects both 
by infall into haloes and through the constant merger of 
gas-containing haloes. A robust power-law mass spectrum 
for the SGMCs is found in our simulations, wherein, in the 
mass range from « 10^ - 10^"Mq, dN/dM oc M^-^^^ ^ 
This is exactly the spectrum that is predicted to account for 
the formation of globular cluster systems around all types 
of galaxies (HP). Further independent support for our con- 
clusions comes from studies of chemical enrichment in hi- 
erarchical formation of the galactic spheroid, in which an 
excellent match to the observed metallicity distributions of 
globular clusters in possible provided that their progenitor 
clouds have a mass spectrum of the same form as our own 
SGMC spectrum (Cote et al. 2000). 

The properties of the individual SGMCs were also de- 
termined in our study. The clouds appear to be in near virial 
equilibrium. Higher resolution simulations will ultimately 
provide better determinations of their properties, as well 
as resolve the substructure of clumps that are expected to 
be the birthsites of stars and star clusters in the early uni- 
verse. Finally, preliminary indications arc that SGMCs with 
the same mass spectrum, and formed by similar processes, 
occur in ACDM models. These results will be studied and 
reported in a future paper. 
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